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FOREWARD 


This report addresses the research which has been performed 
subsequent to that covered in NASA CR-3378: "Case Studies of Clear 

Air Turbulence Using the Diagnostic Richardson Number Tendency 
Formulation". The fundamental emphasis during both these efforts 
was to develop a technique which would be a significant improvement 
over those currently used for flight planning and avoidance of 
clear-air turbulence. Ideally* it should be both quantitative in 
determining its potential intensity and specific in the location of 
regions of relatively high risk. Furthermore* it should not rely 
on specialized data but be functional using the currently available 
raob (rawinsonde observation) system. 

The encouraging results documented in the earlier report were 
considered compelling enough to warrant a closer look into the 
possibilities of a Specific Clear-Air Turbulence Risk (SCATR) index 
approach to the clear-air turbulence (CAT) problem. Unlike the 
aforementioned investigation* which considered the aspects of fuel 
conservation possible by avoiding sustained periods of flight in 
even light to moderate CAT* this study focuses on several cases of 
documented severe CAT. This was because commercial carriers and 
other potential users of such a SCATR index are concerned primarily 
with passenger and crew member safety. 

This research was conducted by the University of Dayton 
Research Institute for the National Aeronautical and Space 
Administration* George C. Marshall Space Flight Center* under the 


technical direction of Mr. Dennis W. Camp and Mrs. Margaret 



Alexander of the Systems Dynamics Laboratory. Part of the effort 
was performed at the Department of Geosciences. Purdue University. 
Portions of the research will be included in the author's Doctoral 
dissertation. The support for this effort was provided by Mr. A. 
Richard Tobiason of the Transport Aircraft. Aeronautical Systems 
Division. Office of Aeronautics and Space Technology. Contract 
No. (NAS8-34687) . Support for the last year of effort was provided 
by Mr. David H. George of the Central Weather Processor project 
of the Program for Regional Observation and Forecasting Systems. 
Environmental Research Laboratories. National Oceanic and 
Atmospheric Administration. 

The author wishes to thank Mr. Dennis Camp for his foresight 
in providing support for such a controversial approach to the CAT 
problem. Further thanks go to the Air Transportation Association, 
in particular Mr. Carl Knable of United Airlines, for their help 
and encouragement. Finally, thanks to Mr. James K. Luers and Mr. 
Patrick Haines of the University of Dayton Research Institute and 
Professors Phillip Smith and Win-Yee Sun of Purdue University for 
providing a sounding board for many of the concepts developed in 
this report. 
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SECTION 1 


INTRODUCTION 

Clear-air Turbulence (CAT) remains a troublesome problem for 
commercial# military and corporate aviation. This is in spite of 
the fact that the problem has been studied for over 25 years. As 
of 1977# the world-wide economic loss to aviation interests due to 
CAT exceeded *25 million annually (Hopkins# 1977). Because of the 
continued increase in air traffic, it is almost certainly much 
greater than that now. Hence# the need for improved operational 
CAT detection and forecasting techniques# particularly for severe 
CAT# continues. 

CAT should actually be considered a generic term in that it 
represents turbulence occuring in statically stable shear layers 
which need not be cloud-free# but may contain nonconvec tive clouds. 
However# it does rule out the turbulence associated with convectively 
unstable planetary boundary layers (whether clear or not) and deep 
convection# such as thunderstorms# as well as mechanical turbulence 
due to wind flowing over the earth's surface. 

The CAT phenomenon# as defined above# is most frequently found 
near the tropopause within jet stream/internal fronts associated with 
active ex tratrop icai cyclones. The relationship between layer 
thickness and vertical wind shear has been statistically related to 
the probable existence of turbulence within a given layer (e. g. # 
Essenwanger and Reiter# 1969). Other observational studies have 
sought CAT in these regions as well as in the lower stratosphere 
<e. g. # Dutton# 1969# Starr and Kemp# 1974# Reed and Hardy# 1972# 
Kennedy and Shapiro# 1975# 80) in order to develop a picture of its 
relationship to both synoptic and meso-scale structures and dynamics. 
Theoretical <e. g. # Lindzen and Rosenthal# 1976# Lalas and Einaudi# 


1 


1980) and empirical evidence <e. g. . Parks et al. » 1984. Hardy et al. , 

1973a) have both pointed to the association of CAT with Kelvin- 
Helmholtz (K-H) instabilities. The spontaneous development of K-H 
instabilities requires the existence within the layer of a critical 
gradient Richardson number> Ri Cf _ . where Ri is defined as 


R 


L 


fo I .?e J 


c ?\v 
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It has also become apparent that atmospheric gravity waves 
must play a supporting role in the development of K-H instabilities 
and their attendant turbulence by locally reducing Ri . The most 
frequent and well-known factor in the excitation of atmospheric 
gravity waves is significant terrain. One particular response to 
terrain has become known as "mountain waves" <Clark and Peltier* 
1977. Smith* 1977). More recently there have been investigations 
into the role of organized convective systems* such as squall lines* 
as barriers to the mean flow* thus acting in much the same way as 
terrain to excite atmospheric gravity waves (e. g. * Keller et al. * 
1983). 

It has become apparent that CAT plays a significant role with 
respect to larger-scale dynamics. It acts as an internal sink of 
kinetic energy when viewed in terms of the large-scale circulation 
of the atmosphere (e. g. * Schenk. 1974* Smith* 1980). While values 
of kinetic energy dissipation rates on a globally averaged basis 
are fairly small (perhaps 2 to 6 W.nf ^ >. they can be much larger 
for individual synoptic scale systems (e. g. * Smith and Adhikary* 
1974; Ward and Smith* 1976; Kung* 1977). Local concentrations of 
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CAT within turbulent shear lagers at jet stream altitudes can 
produce dissipation rates more than an order of magnitude larger 
than the globally averaged values (e. g.» Kennedy and Shapiro* 

1975,80} Dutton, 1969). 

A comprehensive discussion of the turbulent energetics of 
turbulent layers and their relationship to larger— scale structures 
is covered in Dutton (1971). A somewhat less mathematical review 
is contained in Hardy et al. (1973). Aside from its role as a 
dissipative mechanism, however, is the role CAT plays in generating 
potential vorticity. This aspect was studied by a number of workers 
in the early 1960s (e. g. , Staley, 1960). More recent studies (e. g. , 
Gidel and Shapiro, 1979} Holopainen and Nurmi, 1979} Shapiro, 1976, 
1981) have more thoroughly investigated the intriguing implication 
that turbulent processes may actually drive the grid-scale dynamics 
in certain situations. This concept is similar to that advanced by 
Smith (1973) and Smith and Adhikary (1974) using kinetic energy 
budget analyses. 

A discussion of the importance of CAT to aviation interests is 
included in Dutton (1971). Within the economic imperatives of that 
time the concerns were generally in terms of passenger and crew 
member safety and comfort. The advent of the "energy crisis" later 
in that decade added the additional concern of CAT's affect on 
aircraft fuel consumption (e. g. , Winer and Wesler, 1981). In terms 
of flight planning, techniques for location and short-term forecasting 
of CAT for airline operations have remained quite primitive} relying 
on subjective "rules-of-thumb“ and labor intensive meso-scale hand 
analyses (e. g. , Hopkins, 1977} Keller, 1985). 

There seems now to be a strong consensus within the aviation 
community that the CAT forecast currently provided by the NWS is not 
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adequate (Camp et al. > 1980). Despite the great expense involved, 

most airlines which continue to support their own weather staffs 
create CAT forecasts independently of the NWS product. Those 
carriers which no longer have their own weather staff rely mostly 
on pilot reports (PIREPS) to give them some guidance. The decreasing 
volume of PIREPS and increasing number of airlines which are reducing 
or eliminating their weather support point to the need to develop an 
automated, quant i tat ive. and specific CAT risk index which can. at 
least, act as a reliable forecasting tool to aid in flight planning. 
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SECTION 2 


PHYSICAL BASIS 


In this formulation CAT is viewed as a manif estation of internal) 
frictional dissipation of total energy following a parcel of air of 
unit area spanning a well-defined turbulent shear layer of potential 
temperature thickness Ad. The rate of total specific energy 
dissipation within the parcel) vertically integrated through Aq . can 
be shown <Roach) 1970) to be 



where in isentropic coordinates 


<1) 
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( 2 ) 


Avis the magnitude of the horizontal vector wind difference through 

Ae, and an L subscr1pted variable Is defined as: 

-&&/Z 


A more thorough development can be found in the Appendix. The 
development in the Appendix is done using isentropic coordinates 
which yields somewhat cleaner derivations than the pressure 
coordinates used by Roach in his development. The use of 
isentropic coordinates has other benefits for this study; in 
particular) with respect to the analysis of the data to resolve 
the structure of internal fronts. 
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ORIGINAL PAGE ’»2 
OF POOR QUALfTY 


The quantity represents the grid-scale forcin ay larger- 
scale dynamical processes which are attempting to chan the local 
gradient Richardson number, Ri^, within a given layer It is the 
grid-scale component of the so-called diagnostic Rich son number 
tendency (DRT) following the flow (see Appendix), 



no CAT 
CAT 


(A) 


where is the subgrid-scale component. Roach's fun sntal 
assumption is that turbulence occurs within the layer c* a response 
to these forces when they are acting to decrease Ri^ within the 
layer. Turbulence* then* works against these forces :n order to 
maintain Ri^ at a small but nearly constant value. At least one 
recent observational study <e. g. * Kennedy and Shapiro* 1975) has 
shown that the rate of energy dissipation is nearly equal to the 
rate at which the larger-scale deformation field* associated with 
the thr ee-d imens ional meso- and synoptic scale dynamics of the 
jet str eam/fr ontal layers* is acting to increase the vertical shear 
within the layer. Since Ri^ is highly dependent on the vertical 
shear* this suggests that the basic assumption is not unreasonable. 
The equation relating 2^ and can be viewed as a transf ormat ion 

of the grid-scale energy contained within the larger-scale 
deformation field into turbulent energy* rather than acting to 
reduce Rij^. No such relationship is assumed to be relevant to 
layers where these large-scale processes are attempting to increase 


Ri, 


As individual parcels move along through the atmosphere they 
are acted upon by the constantly changing <5 associated with these 
meso- and synoptic scale systems such that Ri^_ , in the absence of 
CAT should increase and decrease according to $ . However, the 
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fundamental assumption suggests that Ri^» uihen averaged over 
and synoptic time scales* would be observed to remain nearly 
within well-established turbulent shear layers provided that 
of sufficient spacial and temporal resolution were available. 


meso- 

constant 

data 
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SECTION 3 


APPLICATION USING ACTUAL DATA 

Bosart and Garcia (1974) used a special analysis data-set 
in an attempt to calculate changes in Ri^ following the flow. Raob 
data were available with sufficient temporal resolution so that fields 
of Ri^_ could be calculated at three-hour intervals. The vertical 
resolution was 2°K (50 - 200 m). Fields of Ri^ were then "predicted” 
three hours ahead using following hand analyzed parcel trajectories 
and compared to the Ri^ fields calculated using the actual analyzed 
data. Many trajectories were calculated for parcels which passed 
through strong baroclinic zones which likely contained layers of- CAT. 

Though they claimed good agreement between the predicted and the 
actual Rij_ fields* their results were severely criticized by Roach 
(1975). Roach contended that the claimed qualitative agreement 
between the observed and calculated Ri^was highly questionable 
first by inspection and second because the contribution of subgrid- 
scale eddies* seen in the complete expression (4) for DRT and not 
included in S' » can not be accounted for in the change of Ri. 

ia 

when parcel trajectories pass through regions of turbulence. The 
accurate calculation of Ri^ following parcels requires very high 
resolution data in both space and time in the first place. When 
turbulence and other subgrid-scale eddies are acting to modify Ri^» 
resulting changes in the detailed vertical structure of horizontal 
wind and potential temperature make defining particular turbulent 
shear layers virtually impossible. 

From the above considerations it becomes clear why calculations 
of Ri^ following the flow are not likely to be useful. Why* then* 
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can it be expected that energy dissipation rates> based on the grid- 
scale forcing of Ri^, be expected to provide useful information on 
the potential for CAT? The reason is simply a difference in the scale 
for which they are appropriately applied. Ri^ is a parameter which is 
highly sensitive to the vertical increment over which it is calculated. 
For a sufficiently detailed data set. calculations of Ri^ reveal the 
tendency for CAT to be correlated to small Ri^ values (< 1). Very high 
resolution data are required in order that the scale associated with CAT. 
for which Ri± has meaning, can be resolved. On the other hand, the 
validity of (j) is not restricted to data having very high resolution. 

V 

Its basis is the role of the grid— scale dynamics in attempting to 
decrease grid-scale Ri- (Ri^ rather than Rij_>. It is reasonably assumed, 
then, that subgrid-scale shear layers, perhaps associated with gravity 
waves, and having smaller Ri (Ri^) will necessarily decrease as well. 
Since the sublayers begin with smaller Ri^ they will reach Ri Cr well 
before the grid-scale layer ever could. 

Indeed. Roach < 1975 ) asserted that the directional variation of 
the horizontal wind vector through turbulent shear layers, associated 
with gravity waves would cause unacceptable errors in the evaluation 
of 2) This is because even ultra-high vertical resolution wind 

data (such as Jimsphere/Jimsondes ) is still not sufficient to resolve 
gravity waves in space and time to accurately calculate their 
contribution to the deformations contained in (JT (eg.. Keller. 1981 ). 
Hence, while it is important that the data used to calculate^be 
analyzed to accurately represent the details of meso-scale frontal 
zones, data fields should be smoothed sufficiently to eliminate 
spurious details. 

The gridded data used for the calculations in this study came 
from raobs. These data were objectively analyzed on either pressure 
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surfaces (as in the case of NMC analysis grids) or isentropic surfaces 
(as in the case of the PROFS/CWP analysis grids) so as to maximize 
the resolution of the meso-scale frontal structure. Although gravity 
waves should be expected to influence the vertical structure of the 
horizontal wind field, particularly in strong vertical shear layers 
where CAT might be expected# analyses of raw raob wind profiles have 
revealed only minor directional variations through even the strongest 
shear layers. None of the quite extreme directional variations 
hypothesized in Brown's (1973) paper# in which he formulated a 
technique to correct for their existence when calculating JT # have 
yet been observed when examining vertical profiles of horizontal 
wind from raw raob data. In addition# these data were smoothed in 
order to remove noise in the wind field that might cause potential 
problems in the calculation of O as discussed above. 

However, energy dissipation rates calculated using data 

linearly fit through a meso-scale layer, will usually result in 

the smoothing or "blurring" effect which increases with increasing 

calculation increment. In his original paper Roach (1970) 

essentially attempted to correlate large values of <5 to pilot 

- Lf 1 

reports of CAT. Roach's case studies proved inconclusive. He 
attributed this largely to the poor vertical resolution of the 
forecast model output data and inadequate analysis data he was 
forced to use. This apparently lead to excessive smoothing of the 
"signal" and inconclusive correlations with CAT encounter reports. 

The data which will be used for this study had significantly greater 
vertical and horizontal resolution than that used by Roach (1970). 
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SECTION 4 


RESULTS 


Estimates of turbulent energy dissipation rates using (1) have 
been calculated for several cases of documented CAT. The documented 
CAT cases which were investigated include: 

(1) Hannibal/ MO accident 3.4 April 1981 

(2) Wisconsin/Illinois outbreak 27.28 November 1984 

(3) Appalacian range outbreak 5.6 February 1985 

The first four sets of figures involve case study (1) which 
used hand-analy zed raob data. This case involved a well publicized 
encounter by a DC-10 with severe turbulence which occured at 
approx imately 0135Z 04 April 1981 near 37 000 ft (225 mb) in the 
Hannibal. MO area. Figure 1 shows horizontal wind vectors and 
isotachs at 250 mb constructed using NMC analysis data. These data 
are available with only 50 mb resolution between 300 and 100 mb. 
While energy dissipation rates can be calculated using this 
relatively coarse resolution, it was of interest to investigate 
the effect of the vertical resolution of the horizontal wind field 
on the calculation of total energy dissipation rates for this case. 
To perform this study, raw raob wind sounding data were hand 
analyzed for a large region of the central U. S. in order to more 
real istical ly represent the vertical shear of the horizontal wind 
for the 50 mb layer centered at 225 mb. 

Figure 2a compares the vertical structure of the horizontal 
wind between 300 and 100 mb at an NMC analysis data grid point 
(broken line) to that for a raob station (solid line; Amarillo. 
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Figure 


Horizontal wind vectors and isotachs <m/s> for 
OOZ 04 April 1981 at 250 mb. 


ORIGINAL PAGE IS 
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Figure 2a Vertical structure of the NMC analysis and raob 

horizontal uiind components for an NMC grid point 
near a raob station for 00Z 04 April 1981. 

Dotted line* labeled^ « shows the direction 
angle of the raob mind vector. 
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Figure, 2 b 


Horizontal ageostrophic vector wind field (NMC) 
for the 250-200 mb layer 00Z 04 April 19B1. 
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TX) virtually at the same point. It can be seen that the analysis 
procedure greatly smoothes the shear in the horizontal winds. Also 
shown is the direction for the raob wind# OC (dotted line). It is 
interesting to note that although the wind amplitude changes greatly 
through this range of altitude# the wind direction changes very 
little. The vertical shear of the horizontal wind is due almost 
entirely to changes in magnitude rather than direction. Inspection 
of numerous other raob wind profiles has shown that this behavior in 
the wind vector is by far the rule rather than the exception. 

The methodology used in this case study was to reconstruct the 

horizontal vector wind field for the 229 mb level using the vertical 

shears for the 50 mb layer between 250 - 200 mb obtained from the 

raw raob profiles. It was the intent of this exercise to more 

approriately represent the vertical shear through the layer while 

preserving the nature of the horizontal wind field. Since the 

structure of the ageostrophic wind field is important in the 

evaluation of (1>« it was felt that preserving this structure was 

evidence that the horizontal wind field had been reconstructed 

properly. Figure 2b shows the ageostrophic vector wind field 

calculated using the NMC winds. The circled 'x' indicates the 

approximate location of the CAT encounter. Figures 3a and 3b show 

for the NMC and raob winds the kinetic energy content# 1/2 V -w 

(heavy solid lines)# and generation# (thin solid lines). 

Energy content rather than specific energy per unit area were 

obtained by multiplying the specific kinetic energy values by Ap/g. 

It is clear that the horizontal structure of the reconstructed 

fields (raob) is very close to those for the NMC winds. Finally. 

figures 4a and 4b show the total energy content dissipation rates 

-2 - 2 . 

(heavy solid lines xlO W*m ) which result using the NMC and the 
raob horizontal winds# respectively. In this case at least# a 
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ORIGINAL PAGE 13 
OF POOR QUALITY 



Kinetic energy content (heavy lines* J. m" 1 ) and 
kinetic energy conversion (thin lines* W» m' 2 - ) 
for the 250-200 mb layer OOZ 04 April 1981 using 
NMC analysis data. 



As 3a except using hand analyzed raob data. 
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Total energy content dissipation rate (heavy 
lines* 10 _i W»nT 2 ) and kinetic energy content for 
the 250-200 mb layer OOZ April 1981 using NMC 
analysis data. The location of the CAT accident* 
at 225 mb* is indicated by the circled 'X'. 



As 4a except using hand analyzed raob data. 


ORIGINAL PAGE IZ 
OF POOR QUALITY 
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more realistic representation of the vertical structure of the 
horizontal wind resulted in a dramatic change in the calculated 
total energy dissipation rate. 

More recent calculations of specific energy rates have been 
performed using finer resolution PROFS/CWP grid point analysis data 
These data come from the standard raob network (North America) but 
have been analyzed to retain much of the detail lost during the NMC 
analysis procedure. These data were also analyzed on isentropic 
surfaces at PROFS/CWP. The horizontal grid is a stereographic 
projection with same geometry as that used by NMC; the resolution 

is about 73 km at mid latitudes. The vertical resolution starts 

o c 

out at 3 K around 700 mb and becomes 9 K near the tropopause 

level. 

Figure 5 shows a measure of the total energy dissipation rate 
(actually* logE^ + 7: Specific CAT Risk or "SCATR" index) calculated 
using OOZ 28 November 1984 analysis data for the 352 K surface. 

While not shown* there were pilot reports of CAT at this level over 
southern Wisconsin and northern Illinois near this time. There was 
also one reported over the state of North Carolina which may be 
associated with the area of large values of £. seen in the lower 
right-hand corner of the figure. The line running diagonally across 
the figure between the two 'x's corresponds to a hypothetical flight 
path between Omaha* NE and Detroit* MI. Figure 6 then shows the 
vertical cross-section of the SCATR index along this path. 

Figure 7 shows the SCATR index field for 12Z 6 February 1985 
at 25 000 feet. Locations of pilot-reported encounters with 
moderate to severe CAT occurred between 20 000 and 30 000 feet 
within two hours after 12Z are also shown. Figure 8 shows the 
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pi^ure. 5. "SCATR index" on the 352 K potential temperature 
surface calculated using PROFS/CWP experimental 
isentropically analyzed raob data for 002 28 
November 1984. The line running diagonally across 
the figure corresponds to a hypothetical flight 
path between Omaha, NE and Detroit MI. 
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Figure 6. Vertical cross-section corresponding to the 
hypothetical flight path shown in figure 5 



Figurt 7. SCATR inda* fiald at 25 000 faat 12Z 06 February 1985 
and modarata-savara CAT ancountars raportad naar thi* 
altituda and tima. Tha diagonal lina raprasants a 
hypothatical flight path batwaan Huntsvilla, AL and 
Hartford, CN. 
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vertical cross-section for the hypothetical flight path indicated 
by the line labeled by Huntsville and Hartford in Figure 7 along 
with the range of altitudes for which CAT was reported as well as 
smooth air reports. 
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6-FEB-1965 122 BASE 10 LOG OF DISSIPATION BY TURBULENCE TIMES 1E7 0 



Figur# 8. Vertical cross-taction corresponding to the hypothetical 
flight path shown in Figure 7. Also shown are the altitude 
ranges and intensities of the reported CAT along this path. 


SECTION 5 


SUMMARY AND DISCUSSION 

In summary, the SCATR index technique has shown some 
promising results. It seems apparent that the improvements 
seen in the case using hand analyzed data, as a result of the 
more realistic representation of the vertical shear of the 
horizontal wind, are also realized in the data analysis used 
in the PROFS/CUP application. Because of the important role 
of meso-scale processes in CAT formation, the performance of 
such a SCATR index should improve as database "enhancements", 
which will increase the resolution in both space and time, are 
brought into even more sophisticated objective analysis schemes. 
These enhancements include profiler and VAS satellite data as 
well as ASDAR and ACAR data. 

The results of the case studies discussed above clearly 
suggest that a SCATR index is not an unrealizable goal. Uses of 
such an index, even in its current prototype level of development, 
are also apparent. These uses include the following: 

(1) "Perfect prog" pseudo-forecasts of SCATR index fields could 

be generated taking advantage of the higher vertical 
resolution aviation forecast winds and temperatures (20 mb) 
soon to become available at NMC for 12 hour forecast 
intervals. This resolution is to be available between 300 
mb and 100 mb. 

(2) Fields of a SCATR index could be made available to NWS and 
other operational meteorologists who then could use it as a 
screening tool to determine which areas might be candidates 
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for applying meso-scale analysis techniques. Both perfect 
prog forecast and analysis SCATR index fields might be useful. 
Forecasters could then qualitatively account for the effects 
of significant terrain (mountain waves)* organized convection* 
and other gravity wave sources 

<3) Another application of such a SCATR index involves its use in 
generating self-consistent flight simulation databases. A 
number of flight scenarios could be created which match CAT 
patterns r ea 1 is t ica 1 ly to the wind and temperatures associated 
with any desired meteorological situation. 

The above applications could be realized within a relatively 
short time* however, there are other applications for which such a 
SCATR index may be useful in the longei — term. These include: 

(4) Use as part of a "MERIT" type short-term forecasting (or 
"nowcasting") model. This could be either by using forecasted 
variables* as for the perfect prog application above* or 
possibly by being incorporated as part of a parameterization 
for internal dissipation. 

(5) The numerical values of a SCATR index provided by a nowcast 
model could then be used as input to flight path optimization 
programs. As mentioned in the introduction* CAT is known to 
produce aerodynamic penalties which increase fuel consumption. 
Passenger and crew member safety would also be improved in 
this way. 

Of course, the best performance of a SCATR index would be 
realized if it is well-calibrated. In order to properly calibrate 
such an index* as well as more thoroughly validate and verify its 
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formulation, it is evident that a field study is necessary. During 
a field study the intensities determined from the grid-scale data 
given by (1) could be compared to the subgrid-scale intensities 
contemporaneously measured by experimental aircraft. 
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SECTION 7 


APPENDIX 


I. The Total Energy Dissipation Rate: Subgrid-scale Information 
from Grid-scale Data 


This development will generally follow that of Roach <1970); 
however* there will be several deviations from his approach. The 
most significant deviation being the use of isentropic rather than 
pressure coordinates. The development will also show more detail 
than what was apparently possible in his paper. 

The equations of motion in isentropic coordinates 


** = f? t <*-7V + H* = -X!f -}k xv + r 

f 9 = c f(fAT 


H = ^ 

3>t 


<A4) 


where - c w T + gz » T = S(p/p 1000 mb* < = R/c^i 


, 1 , -) -r— 

If* -= - )l J ~,\T»jr l turbulent 

' force* 


momentum flux stress 


H = — 


ea r 


turbulent diabatic heating. 


The gradient Richardson number in isentropic coordinates is 



( A5 ) 


The diagnostic Richardson number tendency following the flow is 



(A6) 


The term Dlnd/Dt = 0 in the adiabatic case. For this study* 
where turbulence is considered as a diabatic heat source or sink* 
it will not merely be dropped. The equation of state* p » o RT * 
gives ' 


^ — -Tt. I — (2 


<A7) 


Poisson's equation gives 


At ~ + K-fap p A 
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since l/<fp)Dp/Dt is considered small enough to ignore (small 
compressibility) for the adiabatic case giving Din/? /Dt *=-0. th< 
same assumption will be made for the diabatic case so that 


and 
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This is a form of the Boussinesq approx imation which allows 
density variations only due to diabatic heating. The diagnostic 
Richardson number tendency (from now on called DRT) becomes 


- ft 
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4 %) -JWf;/ 
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( A8) 


What is needed is an expression for DRT in terms of the 
governing primitive equations <A1)-<A4). The turbulent stress 
terms will be retained in order to obtain an explicit form for 
the subgrid-scale DRT as well as the turbulent kinetic energy 
dissipation rate* 

The second right-hand-term can be obtained as follows: 

(Al> gives 

M ) + + ♦ If 
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Using the hydrostatic equation <A2)# and the Poisson equation* 
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Eq (AIO) becomes 
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An expression for the first right— hand— side term in (A8) can 
be obtained directly from the continuity eq. (A3) 
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The complete equation for DRT can be written* for the case when 
when subgrid-scale eddies and turbulence exist* as 



<A13> 


where 
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Cgi7\v 


( 2 > 



(A14) 


is the grid-scale contribution to the diagnostic Richardson 
number. The second term on the rhs of (A14) corresponds to 
shearing deformation by the ageostrophic grid-scale motion field 
which causes tilting of isentropic surfaces. The first and last 
terms, when taken together along the “natural" coordinate of the 
mean flow, represent stretching deformation. Mhen using a non- 
cartesian coordinate system the first term on the rhs must be 
evaluated using the more general form 


( \Y- V ) \V 



*\\/ + 



2 . 


The subgrid-scale contribution is 

£ — 2 . 1 I 
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In the absence of turbulence, then, changes in Ri are entirely due 


The role of the grid-scale dynamic processes in the generation 
of subgrid-scale CAT is considered with respect to its role in the 
dissipation of total specific energy within a parcel as it moves 
with the flow. Consider the layer pictured below having a potential 
temperature thickness A© and horizontal wind speed difference 
through the layer of Av. 
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Figure A-l 


The total grid-scale specific energy of the layer of "depth" A© is 



( A 16 ) 


The figure portrays a well-defined turbulent shear layer specified 
for the limits N* is the Brunt-Vaisala fequency squared* a 

measure of the layer's static stability. In order to evaluate the 
integral in (A16) the structures of v and in the vertical direction 
shall be approximated by a simple linear variation. The author has 
evaluated this integral using a hyperbolic tangent structure for v. 
This exercise yielded virtually the same result as that which will 
be shown using the linear profiles. 


For 
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For let 

t = K + 

and using the hydrostatic equation (A2) gives 
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where Ri^-fccte - Si* /^) 
number* so that f 0 & 6 


local gradient Richardson 


Finally* 



(A16) becomes 
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Now it is necessary to determine the change in following 
the parcel motion at the center of the layer. The depth of the 
layer &B will be kept constant. Although this means that^p* 
and hence the mass of the air volume* will not remain strictly 
constant it will be assumed that it will not alter the basic 
result. The rate of change of e^ can be separated into the 
following components. 



( A2 1 a ) 


(A21b> 


where the subscript 'o' * as before* refers to quantities at the 
center of the layer. 

The first component contributing to De /Dt is 



( A22) 


The right-hand-side of (A22) is zero for steady, nonturbulent 
flow. Of course, this may not be assumed for this situation. 
However, (A22) relates to the changes in e L at the center of the 
layer. It will be assumed that at the center of the layer the 
turbulent heat and momentum fluxes are at a maximum, so that 



The second component of (De^/Dt), involves the transfer of 
energy from the grid-scale deformation field to subgrid-scale 
eddies and turbulence. Equation <A21b> may be rewritten 



and (A13) becomes 
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Equation (All) gives 
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The final expression for (A21) can now be written 


© 


( A25) 


( A26) 


( A27 ) 


37 


( A28) 





(A29) 


can be considered the "explicit" form of the total energy 
dissipation rate. The total specific energy* from <A20)» can be 
written 



(A30) 


Focusing on the last three rhs terms in (A28)> which involve the 
ectual transfer of energy from the deformation field to turbulence* 
in order for Ri^ to remain at a constant (and presumably small) 
value requires a balance between the last two terms such that 


2 . 

A. 




( A31 ) 


If ( A31 ) were not true (at least for established CAT layers)* that 
part not balanced would act to change Ri^. As mentioned in the 
main text* observations (Kennedy and Shappiro* 1975) suggest that 
the balance implied by (A31) may be a reasonable assumption. 


II. The Diagnostic Richardson Number Tendency: Grid-scale and 
Subgrid-scale Effects 


It should be emphasized that (A31) is strictly relevant only 
once the layer pictured by Figure A-l has already become turbulent. 
Consider the general situation. Parcels will pass through regions 
where large-scale processes are acting to increase Ri^ ; and 

hence no turbulence* and then back into regions where tne reverse 
is true: J[ L > O* and turbulence may form if Ri c reaches Ri<* r . The 
general equation governing changes in Ri^ for a parcel is Eq. (A4) 
of the main text or* again* 
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<A32a> 





no CAT 
CAT 


<A32b> 


A separate equation governs the change in Ri L * following the motion 
at the the center of the layer* depending on whether or not CAT has 
formed. Dosart and Garcia (1974) were able* of course* to use only 
<A32a) to calculate changes in Ri^ following the flow. 

The explicit form of the DRT for the turbulent case (A32b) can 
be derived using the explicit forms of £ and]? from (A31) and (A19)< 
respectively. Thus* 4 ^ 



( A33) 


( A34a > 
< A34b ) 


Considering a turbulent sublayer which forms* dissipates and then 
reforms at least several times per hour* when averaging over a meso- 
scale time period (A34b) is the appropriate form. This means that 
for time scales significantly greater than those normally associated 
with the CAT mechanism* it can be seen from <A33> that Ri will be 
maintained at a value of 1/2. 
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